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Abstract: By coupling controllable quantum systems into larger structures 
we introduce the concept of a quantum metamaterial. Conventional metama- 
terials represent one of the most important frontiers in optical design, with 
applications in diverse fields ranging from medicine to aerospace. Up un- 
til now however, metamaterials have themselves been classical structures 
and interact only with the classical properties of light. Here we describe a 
class of dynamic metamaterials, based on the quantum properties of cou- 
pled atom-cavity arrays, which are intrinsically lossless, reconfigurable, and 
operate fundamentally at the quantum level. We show how this new class 
of metamaterial could be used to create a reconfigurable quantum superlens 
possessing a negative index gradient for single photon imaging. With the in- 
herent features of quantum superposition and entanglement of metamaterial 
properties, this new class of dynamic quantum metamaterial, opens a new 
vista for quantum science and technology. 
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1. Introduction 

By offering material properties beyond that which occurs in nature, artificially engineered meta- 
materials are of intense interest. Typically fabricated with periodic features spaced closer than 
the operating wavelength, the system acts as a homogenous material. The earliest introductions 
were in the area of magnetic resonance imaging where conducting elements were used to pro- 
duce artificial magnetism ||TJ [2l [S] (H [51 . Negative index materials (NIMs) with simultaneous 
negative permittivity and permeability have also been engineered. A remarkable property of 
such negative index materials (NIMs) is the possibility of negative refraction (61 . Negative re- 
fraction arises due to the ability to design an interface between concave and convex surfaces in 
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Fig. 1. A reconfigurable quantum superlens built from a cavity-array metamaterial. By 
allowing all-angle negative refraction and evanescent wave enhancement, the superlens 
surpasses the diffraction limitation of conventional lenses. Inset: Electrostatic tuning of the 
intracavity atoms (solid circles) in the cavity lattice provides dynamic control over the light 
guiding and resonance properties of the lens. 

the material bandstructure. This interface leads to a novel dispersion relation not found in nat- 
urally occurring materials, namely all-angle negative refraction (Tl. Simply put, conservation 
of energy forces a reversal of the wave- vector in one dimension, whilst preserving the wave- 
vector in the other dimension. Typically, the appropriate bandstructure engineering requires 
control over the electromagnetic properties on scales significantly smaller than the wavelength 
of the radiation, although any system that realizes the appropriate bandstructure engineering 
will exhibit negative refraction. NIMs with negative refraction have been demonstrated in the 
microwave regime in structures that consists of interlocking metal strips and conducting split- 
ring resonators |[8l[9l[T0l[TTl. By applying transformation optics, metamaterials become a rich 
platform for the control of electromagnetic waves. A striking consequence of such control has 
been the realization of an invisibility cloak |[T2ll . Apart from electromagnetic metamaterials, 
acoustic lEllTll and seismic I15il metamaterials are also areas of intense research. 

Here we introduce a new class of metamaterials distinguished from conventional metamateri- 
als in that it operates in the quantum regime and is easily reconfigurable. The medium, a cavity 
array metamaterial (CAM), comprises of a network of coupled atom-optical cavities. Under 
rotating-wave approximation and in the tight-binding regime, the medium can be treated with 
the Jaynes-Cummings Hubbard (JCH) model |[IS|T7l[l8l[I3[2i|2lll22lEl. Atom-optical 
cavities have already been shown to demonstrate such quantum effects as quantum collapse 
and revival (241, Rabi oscillations 1251 , photon blockade 1261 , and electromagnetically induced 
transparency ll27l . The coupling of these cavities have theoretically been shown to exhibit quan- 
tum phase transitions |[T7]|23l and quantum chaos ll28l . They have also been proposed to act 



as Q-swtiches l|29l , single photon switches [30|, many-body 131] and semiconductor simula- 
tors ||23]| . Importantly, by electrical control of the atomic resonances, the metamaterial prop- 
erties can be dynamically varied, a significant development compared to the more usual static 
implementations. The term quantum metamaterial was first introduced in the context of su- 
perconducting charge qubits inside a superconducting resonator 1321 . Here we extend the con- 
cept of a quantum metamaterial by providing an alternative quantum platform based on cavity 
quantum-electrodynamics. This platform goes beyond classical electromagnetic properties and 
allows the exploitation of a quantum plurality of metamaterial properties. 

As an illustration we present the design for a reconfigurable quantum superlens based on 
the CAM. We envisage a configuration depicted schematically in Fig. [T] where the JCH system 
is manipulated to produce a perfect image using single photons. The aim here is not build a 
better superlens, but to use the superlens as an casestudy of how CAMs can be used to exhibit 
metamaterial properties such as negative refraction. However as with photonic crystal (PhC) 
implementations, CAMs have the advantage of being ideally lossless and do not require the 
operating wavelength to be larger than the constitutent element spacing. Further, because the 
transition energy of each atom in the system can be individually controlled via a Stark shifting 
control voltage, CAMs have the distinct advantage of being highly tunable and reconfigurable: 
features not possible using conventional designs. 

To form a perfect image requires the lossless convergence of the propagating and evanescent 
light components. Conventional lenses only focus the propagating fields, and so the resolution 
of the image is fundamentally limited to features greater than the optical wavelength. Sub wave- 
length features are carried by the high spatial-frequency components encoded by the evanescent 
fields. The loss of the evanescent components lead to the diffraction limit. Near-field scanning 
optical microscopy overcomes this problem by scanning a probe in close proximity to the ob- 
ject, but this is often undesirable for applications such as optical lithography and sensing. It has 
been proposed that a lens built from NIMs can produce perfect far field imaging, exhibiting 
all-angle negative refraction (AANR) and evanescent wave enhancement (EWE) [33 1. Because 
of their ability to overcome the diffraction limit, and the lack of optical axis and curved surfaces 
that AANR affords, such a lens is termed superlens ll33|[34l . 

Negative indexing however is not a prerequisite for superlensing. A different class of meta- 
material from that of NIM, is formed by photonic crystals (PhCs) which uses Bragg scattering. 
PhCs have also been shown to exhibit AANR and EWE ||35][36l|37l, and have the advantage of 
low loss. The disadvantage of PhCs as a superlens is that not all the evanescent components can 
be uniformly amplified. In relation to the PhC as a metamaterial, it is of note that since the size 
and periodicity of the scattering elements in PhCs are on the order of the operating wavelength, 
the medium cannot be considered as homogeneous, which is a necessary condition to identify 
a meaningful permeability and permittivity. 

The theoretical developments of superlenses have also been matched by experimental efforts. 
Superlensing has been demonstrated with microwaves in both NIMs ll38]| and PhCs 1391 l40l . 
as well as other platforms such as silver films in visible light BTI . More recently, near-field 
microscopy using a SiC-based superlens at mid-infrared frequency has successfully imaged 
features smaller than the illumination wavelength B2l . 

We will for the first time discuss cavity arrays as a metamaterial. The inter-cavity hopping 
mechanism of the CAM are markedly distinguished from NIM and Bragg scattering, which 
have been the hallmark of metamaterials until now, and therefore represent another class of 
metamaterial. In particular we will discuss the JCH system as a CAM and investigate it as a 
medium for a quantum-based reconfigurable superlensing device. The tunability of the system 
enables post fabrication dynamical control over the focal point of the lens, and by scanning 
over resonances, importantly allows uniform enhancement of evanescent modes, which is not 



possible with passive PhC-based devices. We further show a gradient negative index lens which 
minimizes reflection from the device. 

2. Jaynes-Cummings-Hubbard Hamiltonian 

We consider a uniform two dimensional (2D) periodic array of coupled optical cavities, embed- 
ded with single two-level atomic systems. Our treatment is implementation independent and 
could be realized in any coupled-cavity array system that realizes a two-dimensional square 
lattice topology. However it is useful to consider a concrete implementation. We specifically 
consider a coupled-cavity array realized in a two-dimensional photonic crystal membrane. In 
such a structure the thickness of the membrane effects confinement in the third dimension, 
whilst the photonic crystal structure defines the cavity array. Such a structure has been consid- 
ered previously in the context of photonic quantum emulators (TtIIJoIO, and one-dimensional 
coupled cavity arrays have been fabricated [43]. The cavity array is a series of quantum oscil- 
lators coupled through the overlapping of the photonic modes of adjacent cavities, such that it 
leads to a tight-binding Hubbard-like model. Each cavity couples to its atom via the Jaynes- 
Cummings (JC) interaction. In terms of the atomic (photonic) raising and lowering operators 
(7+, a~ {al^Ur) at site r, the total JCH Hamiltonian reads (/i = 1),. 

^ = ^ecr/cr~ + a)<2jar + j3(cr/<2^ + aJcr~)- ^ Kalas, (1) 

where Y.{r,s) ^^^^ nearest-neighbor cavities, /c > is the hopping frequency 

(requiring that the ground eigenstate to be symmetric and the first excited eigenstate to be 
anti- symmetric also means that /c > 0), e is the atomic transition energy, (O is the cavity reso- 
nance frequency, j3 is the single-photon Rabi frequency, and the rotating wave approximation 
is assumed. The onsite terms can be diagonalized in a basis of mixed photonic and atomic 
excitations called dressed states or polaritons, |±,n)r = sin0„|g,^)^ + cos0„|^,^ — 1)^, with 
energy = no — A/2 ± y/n^^ + (A/2)^, and mixing angle = ^ arctan[— 2y^j3 /A], where 

A= C0-£. 

Using a Bloch state analysis (see appendix) the band structure in the one-excitation mani- 
fold is given by, 

£± = i(ft, + e-/:)±i^(A-A:)2 + 4j32 . (2) 

For the rotated lattice (depicted in inset of Fig. [2a)), K = 4Kcos{kxd) cos{kyd) and for the unro- 
tated lattice (inset of Fig. [2b)) K = 2K[cos{kxd) + cos{kyd)], where d=\ds — dr\l\fl. Equation 
is an exact solution for a periodic lattice. In our numerical simulations we will verify that 
boundary effects are negligible. Note that Eq. ([T]) describes the underlying connectivity of the 
lattice, and not explicitly its geometry. It is possible for geometrically distinct configurations to 
have the same Hamiltonian, and in the reciprocal lattice the same Bloch analysis would hold, 
resulting in equivalent lensing properties. However, the orientation of the lattice relative to the 
interface is of importance for superlensing, as is explained below. 

In the following sections we will demonstrate how traditional metamaterial techniques can 
be applied to the JCH Hamiltonian, to show that a suitably designed CAM can exhibit the hall- 
marks of a superlens, namely all-angle negative refraction and evanescent wave ehancement. 

3. All-Angle Negative refraction 

A useful tool for analyzing light refraction at an interface is the isoenergy map plotted in k- 
space. Given a dispersion relation E{kx^ky), an isoenergy contour defines the curves over which 
the energy is constant. In this representation the gradient of the energy surface is the vector field 
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Fig. 2. Energy band structure of a JCH lattice with air-lattice isoenergy contours in the first 
Brillouin zone. Bold lines are the isoenergy contours for an operating frequency in free 
space (white) and in the lattice (black). White arrows denote the wavevectors k parallel to 
the phase velocity and colored arrows the group velocity v. At the interface, the ky com- 
ponent is conserved, (a): In the rotated lattice there is all-angle negative refraction (yellow 
arrow denotes the incident group velocity, blue arrow the refraction group velocity), (b): 
In the unrotated lattice, there are no operating frequency propagating modes where ky is 
conserved at the interface, and the photon is reflected (green arrow). 



of group velocities, which points normal to the isoenergy contour as illustrated by the colored 
arrows in Fig.Oa). Using Eq. ©, the group velocity = V^^^^ is expressed as, 



[-It . ^ ^ ]V;-g/2. (3) 

AANR of single photons can take place at the interface between free- space and a JCH lattice. 
Consider the band structure of a rotated lattice shown in Eig. Oa) which is superimposed with 
the lattice and free-space isoenergy contours of matching energy. The contours associated with 
free-space are circles with radius equal to its energy (in natural units). By requiring that the 
surface parallel wave vector ky is conserved at the interface, the group velocities associated 
with these contours determine the refraction angle. In the illustration, an incident photon with 
wavevector ki = (ki^j^kij) and velocity vi will couple to an allowed mode of the lattice, and 
propagate with ^2 = (^2,x7^2j) and V2. The refraction angle is 

Or = arctan(tan^i 3;C0t^2,x) 5 (4) 

where /:2,x is given by Eq. ([T3]) . Since the isoenergy contours of the lattice are convex, we have 
negative refraction (Or < 0), and since the lattice contour is larger than the air contour, this 
occurs for all incident angles. 

To converge the light the isoenergy contour needs to be as circular as possible. This occurs 
at the energy band extrema, and so for a sharp focus it is preferable to work as close to these 
frequencies as possible. 

Refraction is not invariant under lattice rotation. This is because there is a change in the air- 
lattice interface under rotation. Erom the conservation of the ky condition, Eig. Ob) shows that 
the unrotated lattice does not exhibit AANR. When AANR does occur, the propagating modes 
in free space can be brought into a focus to form an image on the other side of the lattice even 
with a planar lattice slab. Such a device therefore satisfies the first criterion of a superlens. In 



contrast, conventional lenses that rely on positive refraction must have curvatures to converge 
light. 

An exact numerical simulation to test the predicted negative refraction (Eqs. © and (O) for 
an air-lattice interface requires knowledge of, and is dependent on, the specific light-cavity cou- 
pling mechanism at the interface of a physical implementation. To demonstrate the underlying 
principles discussed above without recourse to a specific coupling mechanism, we conduct our 
numerical simulation at an interface between two JCH lattices, with no loss in generality. 

We use a segmented lattice as our platform (depicted in Fig.Oa)), where our single photon 
source is initialized in the 'source' region, and there is an identical 'image' region which acts 
as the image plane. The sandwiched 'lens' region negatively refracts the excitation and brings 
it to focus in the image plane. The band structures for the source and lens region are shown 
in Figs. Ob) andOc) respetively. For the chosen operating frequency, indicated by the bold 
isoenergy contours, the system exhibits AANR. 
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Fig. 3. Energy band structure with lattice-lattice isoenergy contours, (a): Lattice configura- 
tion for numerical simulations. Band structures for (b) source and image, and (c) the lens. 
The regions are distinguished by their respective atomic detuning A, and negative refraction 
is predicted at the lens interfaces. The parameters are Ai = 0, j8 = 100/c, A2 = —5.21k. 



Due to the dielectric mismatch of the interface, there must be a finite probability of reflec- 
tion. At the interface, the discrete scattering eigenequation can be used to derive the reflection 
coefficient (see appendix), 

1 -cos(/:i,;, + ^2,x) 

and from conservation of energy the transmission coefficient T = 1 — R. Comparing Eq. ^ 
and Eq. ([5]), there is a trade-off between refraction and reflection, i.e, large negative refraction 
is accompanied by large reflection. This trade-off is illustrated in a comparison plot ( Fig.|4]) of 
the refraction angle and reflection coefficient for different incident angles and varying detuning. 

The propagation of the field in the lattice is governed by the Schrodinger equation | V^(0) = 
^^^^|i//(0)). We consider the case when the source is initialized in an equal superposition of 
atomic and photonic modes. It is instructive to use a directional pulse by specifying an initial 
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Fig. 4. Trade-off between reflection and refraction angle as a function of atomic detuning, 
(a) Reflection coefficient, R, and (b) refraction angle, 0/?, as a function of lens atomic tran- 
sition energy, A2 for different incident angles, 0/ (note that same colored curves correspond 
to same incident angle). Large negative refraction is accompanied by large reflection. The 
model parameters follow Fig. 3. 

state with a normalized Gaussian momentum distribution around k\ = (7r/4,±;r/4) so that it 
is incident on the lens at ±45°, as shown in Fig.[3a). The superposition of the two k modes 
manifests in a coherent interference pattern in the j-direction. 

The lens atomic detuning is set to —5.27k*, which by Eq. (|4]), predicts a refraction angle of 
Or = —25°. Superimposing different time instances, the incident, reflected and refracted pulses 
in Fig. Oa) follow the predicted refraction angle and the trajectory predicted by Eq. ([3]), to 
converge at a location on the image plane. The reflection and transmission coefficients are also 
found to be in good agreement with Eq. ([2Ql) . The incident and reflected polariton (±7r/4, ;r/4) 
coherently interfere near the interface to give an interference pattern along the x-direction. 
Note that there is considerable reflection, so that the population density has been multiplied by 
a factor M in Fig. [3a) and Fig.[5tb) for clearer representation. 

An important property of our system, distinct from the existing PhC superlens implementa- 
tions, is the ability to tune the atomic transition energy, e, after fabrication. Such manipulations 
can be achieved dynamically by, for example, a controlled external electric field via Stark shift. 
This control allows one to tailor the dispersion relation, and hence the light guiding proper- 
ties and focal point of the lens. The effect of changing £ is demonstrated in Fig. Ob), where 
decreasing A2 by 0.73/c shifts the focus 56 sites to the right. 

To show how a point-like source would converge with AANR onto the image plane, we 
specify a source with an initial superposition of Gaussian momentum distribution summed 
over k in Fig. [3c). Taking a snapshot of the propagation at time, t = 300/ K, it shows that as 
the point source propagates into the planar lens, all components are negatively refracted, so that 
an image of the point source is successfully formed on the image plane. As with Fig. [3a) and 
Fig. [3b) there is considerable reflection at the lens' interfaces. 

In our examples so far, the lenses are homogenous and there is an abrupt change at the 
interface. As a result, there is considerable reflection such that the total transmission through 
the lens is less than 25%. As expressed by Eq. ©, reflection increases with the greater the 
change in kx. This can be minimized if we provide an adiabatic spatial change of the atomic 
transition energy within the lens, in effect producing a gradient-index (GRIN) structure. In 
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Fig. 5. Negative refraction of Gaussian polaritonic pulses. The pulse is initialized with a co- 
herent superposition of two momenta ki = (7r/4, ±7r/4). We superimposed different time 
instances t of system evolution with each instance labeled with {Kt,M) where the popula- 
tion is multiplied with M for a clearer presentation, (a): The atomic transition energy in the 
lens is A2 = —5.27k'. (b): A2 = —6k. Predicted trajectories are indicated by the arrows. The 
polariton follows the predicted trajectories for incidence, reflection and refraction. Chang- 
ing £2 changes the focal point, (c): A snapshot at time, t = 300/ K, of the imaging of a 
point- like source by negative refraction (Media 1). (d): The lattice implements a gradient- 
index lens by employing adiabatic variations in A2(x), reducing reflection (Media 2). The 
other parameters follow Fig. [3] 



Fig-Etd) the detuning distribution follows the form, 

r (A2 - Ai ) sin2 ( ) + Ai if < y < w, 

A{x) = Ia2 ifw<x' <W-w, (6) 

[(A2-Ai)cos2(^^i^^+^) + Ai ifW-w<x'<W. 

where is the number of sites from the interface, w is the width of the GRIN region and W is 
the total width of the lens. By fine tuning the GRIN region, the level of reflection can be made 
arbitrarily small, although the physical trade-off is a larger lens. The removal of reflection losses 
is an important development, demonstrating fine control of propagation possible in our system. 

4. Evanescent wave enhancement 

The ability of lenses to resolve images is limited by the wavelength of the light source be- 
cause the high- spatial-frequency modes that describe the sub wavelength features are non- 
propagating and do not reach the image plane. To see this, the dispersion relation in free space, 

kx = 'sj ofi — k^, implies that the modes with ky > CO exponentially decay away from the source 
along the x-axis. Existing superlens proposals overcome this diffraction limit by amplifying the 
evanescent wave (EW) components. 



PhC-based evanescent wave enhancement (EWE) devices can be regarded as a type of res- 
onator. At resonance, the transmission of the evanescent components is divergent. The total 
transmission across a lens of width W, derived from taking the summation of the multiple scat- 
tering events at the left and right interfaces is, 



— (7) 

cxp{-2ik2^,W)-Rl^ ' 
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where Tij (Rij) is the transmission (reflection) amplitude at the interface between region / and j. 
At the resonance condition exp(— 2//:2,x^) ~ ^23 ~ ^' transmission is divergent. The resonant 
condition is just the condition for total internal reflection where the accumulated phase shift in 
a round trip is a multiple of 2;r. 

The resonant bound modes 136] allow a build-up of these bound states to produce an ampli- 
fied evanescent tail on the image side of the lens. These resonant modes are discrete, therefore 
they will not amplify all evanescent modes. NIM-based superlenses do not have this limitation, 
but because their fabrication is based on conducting elements they suffer the problem of loss. 

Since JCH-based superlenses can dynamically shift the resonant points, they can overcome 
the limitation of PhC-based superlenses and amplify a contiguous range of evanescent modes 
(but not simultaneously). 

Resonant bound modes can either exist along the interface (surface bound modes) or in the 
bulk (bulk bound modes). As both mechanisms follow the same underlying resonance principle, 
we demonstrate only the latter. 

A quantum equivalent of an EW is the evanescent tail of a stationary state. We prepare our 
system such that the detuning in regions 1 and 3, Ai and A3 respectively, are sufficiently differ- 
ent from the detuning in region 2, A2 (see Fig. Oa)), so that we can setup an eigenstate where 
evanescent tails exists in region 1 and 3. This setup is analogous to that of a square well. Region 
4 is the lens, and region 5 our image plane. 

The resonant condition for EWE can be achieved by tuning the atomic transition energy of 
the atoms in the lens region or the width of the lens. We will only consider the former here, as 
it is more appropriate to the idea of a reconfigurable device. We solve for the time evolution of 
the system and observe in Fig. Ob) that significant coupling to the lens occurs at distinct values 
of the atomic detuning, A4, where the resonance condition is met. 

When the lens is exactly on resonance, as demonstrated in Fig.^c) where A2 = A4 = 0, the 
population exchange occurs between the source and the lens is that of coupled homogenous 
resonators, i.e., P4{t) = sin^(^lr), where is the characteristic mutual coupling. Following the 
increase in P4, the EW incident on the lens is transmitted amplified by 6 orders of magnitude 
after a time ^ = lO^/c"^. 

When the lens is only quasi-resonant (Fig.^d)), the population exchange between the cou- 
pled resonators (source and lens) can be approximated by. 
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(8) 



where T] is the difference in the eigenenergy of the source and lens. By fitting P^it) to nu- 
merical results, we find that for A4 = 0.305k-, t] - 10"^ /c. This is in good agreement with the 
minimum energy difference between source and lens obtained by solving the lens Hamiltonian 
directly. 

Fig. Oc) shows that at t = 10^/c~^ the incident EW is amplified by a third of the exact 
resonance case. Thus, although exact resonance is an optimal condition for EWE it is not a 
necessary condition for enhancement. When the lens is tuned away from resonance, the degree 
of EWE can quickly diminish as seen in FigOe). 
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Fig. 6. EWE by tuning atomic transition energy, (a): Lattice schematics for demonstrating 
EWE. The evanescent wave (EW) is initiahzed in the region 3 by preparing an eigenstate 
of energy — CO ^ lOO/c in the region 1-3. The lens serves to amplify the field via resonant 
coupling, (b): Population buildup in the lens (P4) occurs at the quasi-resonances (peaks). 
(c)-(e): Population profile taken along the x-axis and time evolution of population P4 for 
selected A4 values. Dashed blue line shows the rapid drop-off of the EW without the lens, 
such that the incident population or probability would be 10~^^. The lens provides the 
enhancement on resonance seen in (c) and (d). The parameters are Ai 3 5 = 0, A2 = and 
p = 100k-. 



The diffraction limit restricts the resolution of conventional lenses to the operating wave- 
length, Aq. Our lens' resolution, 5 = 2n/kmax, is determined by the maximum k that still satis- 
fies the resonant condition. However we would also like to resolve all the ^-components leading 
up to kmax- This implies minimizing the bulk energy band spectrum so that the deviation from 
the resonant energy is always small. The drawback of this is a reduction in sharpness of focus. 
A better solution is to introduce surface mode resonance. This can be achieved by having a dif- 
ferent £ at the lens surfaces from that of the bulk. As shown in Fig. [71 the flatter surface mode 
band (see appendix) provides the necessary minimal deviation from resonance to maximize 
^max, leaving the bulk mode to provide the AANR and focal sharpness. 

Fig. |7] shows that ^max = 27t/d — (Oq, because beyond this the evanescent modes fold back into 
the light cone and the associated bound modes become leaky states ll36ll . Thus the maximum 
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Fig. 7. Surface and bulk energy band structures for evanescent wave enhancement. There 
are regions in the bulk energy spectrum which deviate from the frequency, OX). In these gaps, 
EWE via bound bulk mode resonance can not occur. By having surface mode resonance, 
these gaps can be filled. The thin surface energy spectrum around the operating energy, ft^, 
maximizes the range of k for which EWE can occur, whilst the broad bulk energy spectrum 
provides AANR and focal sharpness. Surface modes can be achieved by independently 
tuning the lens surfaces from the bulk (as depicted in the inset). 



resolution of our lens is, 

I— a / Ao 

so that for sufficiently small inter-cavity spacing, d < Ao/2, the resolution exceeds that of 
conventional lenses. 

Inter-cavity spacing is however limited by the size of the cavity. This means that to beat 
the diffraction limit, one at the very least needs the cavity size to be less than Aq. Typically 
cavity resonant wavelength is twice the cavity size, so the subwavelength resolution condition 
becomes, (O > (Oq. The non-linear interaction introduced by the cavity atom allows, beyond 
that which is available through just inter-cavity hopping, the cavity resonance frequency to be 
greater than the operating frequency. 

Using Eq. Q where = (Oo, the relative resolution of our lens can be approximated by, 

(10) 

where we have assumed small detuning and that the minimum possible spacing between sites 
is half the resonant wavelength. Equation ([TOl) gives the factor by which our lens beats the 
diffraction limit. Since 5o/5+ < 1, only the resolution from the negative energy branch, 5~, 
can better the diffraction limit. Conventionally j3 /o), restricted by the so-called fine structure 
constant limit, is of the order 0.01 (although larger values are possible for unconventional cou- 
pling mechanisms ll44l ). so the improvement over the diffraction limit is typically small. 



5. Experimental feasibility and outlook 

The lead time from theoretical inception to experimental realization can be long. NIMs were 
originally proposed in 1968 |6|, but it was not physically realised until 2000 |8|. The case 
for the experimental feasibility of CAMs is fortunately alot clearer. Optical microcavities have 
been created with whispering gallery modes (WGM) B31l . fabry-perot microcavities B6ll47l . 
photonic bandgap (PBG) defects 1481 , and in slot waveguides {49\ . The coupling of PBG 
nanocavities into a one-dimensional array has also been achieved |43l . Recently the coupling 
of Nitrogen- vacancy (NV) centers to PBG cavities ISOl and WGM microdisks ISTl have been 
experimentally verified. With the astonishing advancement in microcavities and arrays of mi- 
crocavities, and their coupling to multi-level quantum systems, it is feasible that the next ad- 
vancement in experimental development would be that of arrays of microcavities coupled to 
multi-level quantum systems, in other words, CAMs. The dynamic control of Stark tunable 
solid-state emitters have been been demonstrated with diamond colour centres 1521, quantum 
dots |[53]| single molecules ll54l . and single-ion dopants ll55]| . 

The development of microcavities can be characterized by their size and quality factor, Q. 
WGM microcavities have experimentally produced 2 10^ in microspheres 1561 and 10^ in 
microtoroids (571 [58l. However their sizes are relatively large: 10^ /im^ and 180 11100? respec- 
tively. PBG microcavities have achieved 2^10^ with cavity mode volume, V ~ (A/2)^, where 
A is the operating wavelength 1591 . NV couplings in PBG microcavities have been calculated 
as j8 ^ 10^^ Hz 1601 . Assuming that photon hopping limits Q, we can approximate the inter- 
cavity tunneling frequency sls K = CO / Q. For the superlensing properties presented in this work, 
j3 = lOO/c, requiring in the visible light regime 2^10^, which is at current experimental limits. 
In PBG arrays with over 100 microcavities however, only 2 ~ 10^ as yet been experimentally 
verified (431. Coupling sub wavelength- sized slot- waveguide cavities in 2D lattice has also been 
discussed (611 . 

Atom-optical microcavities have demonstrated many quantum effects ranging from quan- 
tum collapse and revival to eletromagnetically induced transparency. It is upon these quantum 
effects that CAMs' true potential lies. In particular, a transmission line formed by coupling 
superconducting charge qubits prepared in a coherent superposition of quantum states has been 
studied to show an oscillating band gap and transmittance (32l . First experimental steps to real- 
ize such materials have been reported (62l[63l[64l- Along similar lines for future investigation, 
coupling of a three level A- system (e.g. using quantum dots (65l or diamond defects (66l ). to 
the cavity modes could yield a superposition of dispersion relations and thus create a metama- 
terials with a quantum superposition of material properties. One of the ground states in the 
individual cavity A-system would be strongly coupled to the cavity modes while the other (|/)) 
is not (Fig. [8]). The atom-photon coupling offers intriguing potential for entirely new quantum 
devices. For example, if some of the atoms are prepared in a Greenberger-Horne-Zeilinger-like 
(GHZ-like) state \gg...g) + I//--./), then the system will exhibit a superposition of dispersion 
relations. This in turn implies a superlens with two focal points in quantum superposition. 

6. Conclusion 

In summary, we have combined the previously unrelated fields of quantum mechanics and 
metamaterials, by proposing cavity arrays as a new class of dynamic metamaterial. Operating 
at the quantum level, it opens up new possibilities for quantum optical devices. By applying 
traditional metamaterial techniques we showed that the CAM can exhibit the features of a 
superlens. In a more general sense, this work lays down the framework for local manipulation of 
photons, the quantum superposition of metamaterial properties, the preservation and interaction 
with entangled fields, and other non-local effects, in cavity array metamaterials, creating a 
new area of investigation in quantum transformation optical phenomena. This invites quantum 



Fig. 8. Quantum superposition of material properties with cavity A-system. Each cavity is 
a A-system (e.g. quantum dots or NV centers in diamond) consisting of two atomic ground 
states and one excited state, \e). One of the ground states, is coupled to the cavity mode 
via coupling parameter /3 and the other ground state, |/), is only weakly coupled. When 
some of the atoms are in a GHZ-like state, the metamaterial will exhibit a superposition of 
dispersion relations. The superlens then would exhibit a superposition of two focus points, 
represented by the blue and red arrows. 

technology into the realm of metamaterials. 
Appendix 

1 . The method described here for the derivation of the band structure solution in Eq. (O follows 
Ref. ||23]| . In the presence of intercavity coupling, the onsite energies Ey are no longer the 
polaritonic energies and in general, satisfy the relation, 

Y,Hrs\^,s)=Er\(i>,r), (11) 

s 

where Hys is the Hamiltonian that relates site r to site s. Employing Bloch's theorem for periodic 
structures, 

\(\),s) = |0,r)exp[/^-(J,-4-)], (12) 

where dy is the displacement to site r, and k= {kx^ky) is the wavevector associated with the 
crystal momentum, Eq. (fTTl) becomes an energy eigenequation whose eigenvalues are the en- 
ergy band structure or the dispersion relation of the medium. 

2. For an incident field of energy Ei and wavevector ki = {ki^x^hj) with transmitted field of 
E2 and k2 = {k2,x^k2j) at the interface, energy conservation (£^1 = E2) and phase matching 
(kij = k2j = ky) requires that ^2,x satisfies the condition, 

K=0}2-Ei^l3i/{Ei-£2) , (13) 
where K = 4k*cos(/:2,x<^) cos{kijd) for the rotated lattice. 



3. For the derivation of the reflection coefficient (Eq. ©), the state vector can be expanded in 
the bare atom-photon basis, 

\V) =T,{Cr\8A)r^dr\e,0)r)(S)\8.0)s (14) 

r s^r 

where \g^n)r and \e^n)r denote the ground and excited state respectively, with n photonic 
excitations at site r. Given the symmetry, we consider a 5 -site unit ceU in X-configuration 
that is translational-invariant along the j-direction. Using the standard eigenenergy equation 
J^\\l/) = E\\i/), one arrives at a discrete scattering equation for each region of the lattice and 
the interfaces. In particular in region j with associated parameters (cOj, Kj), we have 

where a given site at coordinate (p^q) is surrounded by four nearest-neighboring sites at coor- 
dinates (/? ± 1 , ^ ± 1) and the conservation of energy requires Ej = E. At the interface centered 
at the origin (0,0), 

k:i(co-i+ci,o) + /C2(co,i+c_i,o) = 0)2 + -r^^ ^ co,o • (16) 

We make the typical assumption that region 1 consists of an incident and a reflected wave 
component, 

where kj^p — kj^x — ky and kj^q = kj^x + Z:^. In region 2, the transmitted wave has the form, 

Cpq = te'^^^pPe'^^^^'^ , (18) 



Cp q = e'^'^pPe'^'^^"^ + re-'^i.pPe-^K,^ , (17) 



where r and t are used here to denote reflection and transmission amplitudes respectively. Sub- 
stituting these solutions in to the interface equation (Eq. ([T6l) ) and applying continuity condition 
^ = 1 + r, we arrive at the reflection coefficient R= |rp, 

Kl^Kl-2KiK2CO^{ki^x-k2,x) 

K= — 7i 7i -. \^^) 

Kf + K-^ - 2K-1 K-2 C0S(^1,;, + k2^x) 

Assuming uniform coupling Kj = K, we retrieve the required expression. Finally since a polari- 
tonic pulse has a momentum distribution G(k), we define an effective reflection coefficient, 

^eff = G{k)R{k)dk. (20) 

J -It J-K 

4. To calculate the surface mode energy band, we need to take two adjacent sites along the 
surface as the primitive cell. After applying B loch's theorem the Hamiltonian is, 

( O) J3 -K[\^Qy.^{iky^/2d)] 0\ 

jS 

-K[\^Qx^{-iky\/2d)] CO P 

V p eJ 

where £s and are the surface and bulk atomic transition energies respectively. We solve for 
the eigenvalues to get the surface mode energy bands. 
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